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We present a general method for the high-teinperature expansion of the self-energy of interacting 
particles. Though the method is valid for fermions and bosons, we illnstrate it for spin one half 
fermions interacting via a zero range potential, in the Bose Einstein Condensate - Bardeen Cooper 
Sclirieffer (BEC-BCS) crossover. The small parameter of the expansion is the fugacity z. Our 
results include terms of order z and z^, which take into account respectively two and three body 
correlations. We give results for the liigh temperature expansion of Tan’s contact at order z’ in the 
whole BEC-BCS crossover. We apply oin method to calculate the spectral function at the unitary 
hmit. We find new structures which were overlooked by previotis approaches, which included only 
two body correlations. This shows that including three-hodv correlations can play an important 


role in the structures of the spectral function. 
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I. INTRODUCTION 

It is a general trend of the physics of ultra cold atoms 
to become the field of experiments simulating condensed 
matter problems [J. In particular, interacting ultra cold 
fermionic atoms can be used to simulate interacting elec¬ 
trons in solid state physics. What is remarkable i.s that 
there is essentially no unknown relevant microscopic pa¬ 
rameters and that some (e.^. the strength of the inter¬ 
action) can be changed at will experimentally. Therefore 
the compfuison to theory can be very accurate. An ex¬ 
ample is the measurement of the equation of state for 
ultra cold ^Li atoms |21[3]. In particular, the high tem- 
pcratiu’c equation of state can be very wcU fitted iLsing 
a virial expansion [2]. At unitarity (where the scattering 
length diverges), the third order virial coefficient deter¬ 
mined theoretically [^EllE] fits perfectly experimental re¬ 
sults. More recently, the method of Rcf.|^ was extended 
to the expansion of the self-energy in the context of the 
2D BEC-BCS crossover [TO]. There the que.stion of a 
depletion (the so-called pseudogap) in the spectral func¬ 
tion, which physically represents the density of states of 
excitation energies after the creation or annihilation of a 
particle of wave vector k was asked. Doing a lowest or¬ 
der calculation, which is a high temperature expansion of 
the T-matrix approximation, the authors of HE!] found 
that such a depletion can be found. However, the pecu¬ 
liarity of the 2D BEC-BCS crossover is that there exists 
always a two-body bound state called a dimer. There¬ 
fore it is quite natiu'al (see Appendix O for a physical 
explanation of this fact) to find a pseudogap feature in 
the whole crossover. This is of course not the case in the 
3£) BEC-BCS crossover at the unitary limit, where there 
is no dimer state in the two body problem. A high tem¬ 
perature expansion of the one-particle Green’s function 
in a trap including 2-body correlations was also studied 
in m- Using the T-matrix approximation, it was found 
at the unitary limit that there exists a depletion of the 


spectral function for moderate values of the wave vector 
k 1121 1131 . Using a different theoretical approach, the 
spectral function was also studied in |14l . Experimen¬ 
tally, the observation of a pseudo gap hiia been foimd 
in R.efs. [T^ llti) . This result was also found in some nu¬ 
merical calculations in Ref.|17| (which are numerically 
difficult since one hits to find a real frequency quantity 
from imaginary time calculations). We recover this be¬ 
havior in section [IIIl because mu' expansion at this order 
coincides with the high temperature behavior of the T- 
matrix approximation. Our goal Ls to go beyond this 
approximation, in a high temperatme controlled calcu¬ 
lation. Indeed most of the analytical calculations of the 
spectral function to the best of our knowledge, have been 
done using this approximation. So, we present in section 
[V| oiu analytical expressions and in section |VI| our mi- 
luerical results beyond this approximation (i.e including 
three body correlations) at the unitary limit. 

The paper is organized as follows. In section the 
general formali.sm is introduced. Then we present results 
of the high-temperature expansion for the self-energy, 
spectral fimetion and Tan's contact to the lowe.st order 
in section 11111 Section II VI is devoted to a virial like cal¬ 
culation of Tan's contact at third order. In section ^ 
we give expressions for the fermion self-energy at next 
order in the fugacity, including 3-body correlations. The 
technical details leading to the analytical expressions are 
explained in Appendices [U] and [H] In sectioti [Vl| we show 
our numerical results for the spectral function including 
three body correlations, pointing out the differences three 
body correlations make compared to the usual T-matrbc 
approach (see Ref. ll2| for instance). Finally, we conclude 
in section rVIIl 

In this work, we will take the Boltzmann constant fcs 
and the Planck’s constant h equal to unity. 
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II. GENERAL FORMALISM 


III. CALCULATION OF 


The general idea of the high temperature expansion 
[5] is to expiind the fermionic non-interacting Green’s 
functions in powers of the fugacity 2 = Here, 13 = 
IjT [T the temperature) and /t is the (Jieinical potential. 
If we denote by r) the non-interacting fermionic 

Green's function, we have 

G™(k,T) = ( 1 ) 

with 0(x) being the Heaviside function, 0 ^( 1 ) = 
l/(e^* + 1) is the Fcrmi-Dirac dLstribution and = 
k^/{2m) (m the mass of an atom), r is the imagi¬ 
nary time. Next, the Fermi-Dirac distribution is ex¬ 
panded in powers of the fugacity : np’(ek ~ 1^) = 

We can write 


P-ki 



FIG. 1: Lowest order diagram for the self-energy. 


G‘°’(k,T) 


n>0 


We have defined 


The lowest order diagram for the self-energy is shown 
in Figj^ Previous authors [Slini] have used this diagram 
in the context of the 2D BEC-BCS crossover. The ana¬ 
lytical expression is 


G<"'"*(k,r) = -0(T)e-"^^ (3) 

G<°''‘’(k,T) = n> 1 (4) 







(5) 


Therefore is retarded, while for n > 1, is 

not retarded. Notice that G*'’ '** does not depend on the 
chemical potential Diagrammatically, since G**^'°* is a 
retarded function, we represent it as a line with tm iurow 
going from left to right if increasing time gots to the right 
(this is the Green’s funetion of a particle in vacuum). On 
the other hand, , for n> 1 . is not retarded, and we 

represent it as a n-times slashed line, which can be ori¬ 
ented from left to right or vice versa. This formalism was 
used in [S| in order to calculate the fermion oecupation 
number n^, whieh is obtained by closing in time the ex¬ 
act Green’s function. By integration on wave-vectors, we 
get the density n(^, T) and by the Gibbs-Duhem relation, 
we can calculate the equation of state P(n, T) where P is 
the pressure of the system. This low fugacity expansion 
is a high-temperature low density expansion, as can be 
the low density expansion in a Lee-Huang-Yiuig like cal¬ 
culation for a moleculiu' BEC |18II11J| . However, here the 
density is low compared to where = y/2K/(mT) 
is the thermal de Broglie wavelength. Instead in [TS) [T5], 
the density is small compared to a~^ where a Ls the scat¬ 
tering length. 

In this work, we use the high temperature expansion 
ill order to calculate the retarded self-energy Ei?(fc,w) 
ill powers of the fugacity z. This allows ils to calculate 
the retarded Green’s function Gi?(A;, w) and the spectral 
function A(A:,tj) = —l/ 7 rG(Gi}(fc,w)). 

This approach was used in and HOI at lowest order 
in the fugacity 2 in the context of the 2D BEC-BCS 
crossover. Our present calculation is instead up to order 
2 ^ for the "iD BEC-BCS crossover. In particular, we take 
into account the 3-body problem. 


T 2 denotes the 2—body T—matrix in the case of a contact 
interaction. In terms of Feynman diagrams, it is given 
by the sum of ladder diagrams [S]. We use the Galilean 
invariance property 72 (F’.r) = r), and we 

denote f 2 (s) the Laplace transform of 72 (0, r). In 3D for 
equal masses, we have t 2 {s) = ^ ~ y'—m.s] . We 

find 


T2iP,T) 



( 6 ) 


where C-y is a Bromwich contour I2UI . 7 is such that 
the integrand is analytic for s such that 3?(.s) < 7 . For 
a~^ < 0 , it is sufficient to have 7 < 0 , while for > 0 , 
7 < —El due to the molecular pole —Ei = —l/(m<t^) of 
t 2 . In Eq.(M, we deform the contour along the real axis 
and we fincT 


T2iP,T) 




+ 


f+ao 

I dx I 
Jo 


V 2 (x) 


(7) 


where 0(x) is the Heaviside step function, Z,,, = 
87 r/(jn^a) is the molecular residue and P 2 {^) = 

—l/ 7 rS(f 2 (i + iO'*')) = ‘l/m^^^y/x/{x + Ei) is the spec¬ 
tral density of the 2-particle 7-matrix. The physical in¬ 
terpretation of the two terms of Eq.Q is clear: the first 
term comes from the two-body bound state of energy 
—El (the dimer), while the .second is due to the two- 
partielc scattering states continuum of energy x. Us¬ 
ing this expression in Eq.(j^. we can easily perform the 
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Foiirior transforru with w„ = 7rT(27i + 1) a This is the result of Eq.(13) of Rcf.[Tn] which was found 

fermionic Alatsubara frequency. We find for a 2D .system. It was aLso derived in Ref. [ 21 ] iu the 

2 eontext of the physics of niagnons. 

iUn) = Z Fi{k,iuJn +^l) + Z Hi(k,iUn + f-t) (8) 



where the two functions Fi and Hi are defined through 




'( 277)3 

e{a-^] 


+ CC 

dx 


P2(x) 




' 4m 


r' /P'^ 17 (P-k) 


am 

( 9 ) 



for a general complex energy E outside the real axis. 
TIutsc expressions call for a simple physical interpreta¬ 
tion. They contain energy denominators of the form 
E — AE where AE = P'^/{4m) -I -1 — (P — k)^/(2m) 
is clearly the energy of an intermediate state with a 
spin up and a spin down interacting particles of cen¬ 
ter of mass kinetic energy P^/(4m), total momcntiun 
P and relative motion energy x, pliLS a spin down hole 
energy —(P — k)3/(2r7i). This can be seen directly on 
the diagram of Figj^ where the two interacting particles 
are represented in the symbol which propagates for¬ 
ward in time, while the hole, represented by the slashed 
line propagates backward in time. Terms of the type 
l/(E-(P'^/{4m)-Eb - (P -k)2/(2m))) of course rep¬ 
resent intermediate states where the relative motion is in 
the two-body boimd state of energy 

From the expre.ssions Eqs.(j^ and ( |10] . we can easily 
perform the analytical continuation to the real axis so as 
to get the retarded self-energy. Indexid, we see that F\ 
and Hi arc analytic functions of E in the upper and lower 
half complex plane. Therefore we find for the retarded 
self-energy 

= zFi{k,uj + ft+ i0'^) + Hi{k,uj + iJ + iO'^) 

( 11 ) 

From this, we come to the conclusion that the diagram 
of Fig{^ actually contains terms of order z and z^. We 
can actually simplify the expression for Fi- We use the 
relation for the function t 2 {z) 

( 12 ) 

Z+Eb Jo (2-.t) 


A. Examples of A(k.ui) 





(m+fi)/Ep 

(c) 


which comes from the analytic properties of t 2 {z) and 
which is therefore valid in any dimension (with the proper 
Then we see from Eq.(|^ that 




4m 2m ^ 
(13) 


FIG. 2: Spectral functions ,4(jfc,ui) versus w-|- ^ at T = 
Tf- (a) l/(fcfo) = 1 (BEC limit), (b) 1 /( 1 : 770 ) = 0 (Unitary 
Limit) and (c) l/(fefa) = —1 (BCS limit). In each panel, 
three values of the wave vectors are shown: k/kp = 0 (black 
solid line), k/kp = 1 (red dashed hne) and k/kp = 2 (blue 
dash dotted line). 
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We compute the spectral function yl(fc, w) from the 
identity 


A(fc,cj) = —G(GR(fc,w)) 

TT 


(14) 


The retarded Green's function is deduced from the com¬ 
puted retarded self-energy llfi{k,ui) via the Dyson equa¬ 
tion 




uj + fi + i0+ - k‘^/{2m) - T,R{k,u) 


(15) 


Figs 2(a)|2(b) and |2(c)| typical examples of 
:tions A(i',tj) using the lowest order self en- 


Wc show in 

spectral functions A(A;, cj) u sing 
ergy of order z in Eq.( |llt. This is the 3D version of 
what wiis done in 2D in I91ll()|. However, contrary to the 
2D case, in 3D there is no moleculiir state for a~^ < 0. 
We define the Fermi momentimi kp in terms of the to¬ 
tal particle density n = kp/{37r^). The Fermi energy 
Sp = kp/{2m) = Tp. It is explained in Appendix [A| how 
to determine the fugacity z for given T/Tp and 1 /[kp a). 

Starting from the BEC side (l/(A:j7a) = 1,2 = 0.17), 
we find that the spectral function has in general two 

a peak on the positive 
u 


distinct features (see Fig 2(a) 
fi part, centered around ui 


fi = A:^/(2m) (essen¬ 
tially a Lorentzian shape) and a broad structme starting 
below —Eb + k^ / {2m). These .structures have a very sim¬ 
ple physical interpretation. The peak centered around 
k^/{2m) corresponds to the creation of a spin f fermion 
of momentum k in a gas of classical dimers. In the BEC 
limit a is small and the scattering of the created fermion 
on the dimers is weak. These enable to find the position 
of the peak, with a mean-field shift and the broadening 
of the peak which can be calculated using a simple Fermi 
Golden Rule iirgrimcnt. The broad structure corresponds 
to the annihilation of a spin t fermion in a classical gas 
of dimers. One is left with a single spin4. fermion in a gas 
of classical dimers. The kinetic energy of this fermion is 
(P — k)^/(2m) where P is the momentum of the dimer. 
This explains the shape of the curve as well as the exis¬ 
tence of a threshold. All the necessary details are given 
in Appendix 

On the BCS side {l/{kpa) = —1, z = 0.7), the struc¬ 
ture is simpler: a peak located around the kinetic en¬ 
ergy k'^/ (2m) (also with a Lorentzian shape) with a finite 
broadening (see Fig 2(c)[ l. Again, this hiis a simple inter¬ 
pretation; the created spin t fermion is going to scatter 
weakly on the spin 4- fermions. As a consequence, there 
will be a mean field shift at the position of the peak. The 
width will be given by a Fermi Golden Rule calculation. 
More detailed arguments can be forind in Appendix [B} 

In the unitary limit {l/{kpa) = 0, 2 = 0.5), a double¬ 
peak structure emerges at small momentum imd it grad¬ 
ually evolves into one peak at large momentum. A small 
bump still exists near the main peak at the Fermi mo¬ 
mentum, which makes it look very asymmetric. These 
characteristics cannot be interpreted simply in the same 
picture as on the BEC or BCS side. Coming from the 
BEC side, the bound state will disappear. Thus the 


corresponding left feature should disappear in the uni¬ 
tary limit and there should exist only one peak on the 
BCS side. Mathematically the double peak structure re¬ 
sults from the real piirt of the self-cncrg^' since the sec¬ 
ond term of the imaginiuy part in Eq.i^ vanishes and 
the remaining term is fairly smooth in the \mitary limit. 
So we attribute this structure to the strong scattering 
between particles in the unitary limit and three-body 
physics or more-body physics may play a significant role 
in the strueture of the spectral function. 


B. High negative frequency, high wave vector 
behavior and Tan's Contact 


For finite ui and k, the dominant contribution is the 
first term in Eq.illi. It is of order z. However, if 


one is interested in evaluating the Tan’s contaet, one 
should consider the high k limit. In this hmit, as noted 
by the authors of Refs. [221 I23|, the dominant contri¬ 
bution to the spectral function or the imaginary part 
of the retarded self-energy comes from the region w ~ 
—k^/{2m). However, in this limit, this is the imaginiiry 
part of the second term in Eq.( |ll[ l which is domintuit. 
Indeed, one finds eiisily, u.sing Eqs.® imd 1(1 1 that 
^[Fi(k,u + ti + i0+)] = + fi + i0+)]. 

In the limit where /Su —> —oo and |lj| '3> |fr|, wc 
find that the term of order z^ is dominant. Therefore, 
the occupied spectral function A_(A;,uj) = f{ui)A{k,ui) 
(/(w) = + 1)“^ is the Fermi-Dirac distribution) in 

this high k limit is dominated by the term of order z^. 
This is shown in Figs 3(a)|3(b)| and 3{c) This hiis as a 
consequence that the Tan s contact is given by the ap¬ 
parently non-dominant term. In order to calculate the 
fermion occupation number n^. one needs to integrate 
on the frequency the spectral function times the Fcrmi- 
Dirac di.stribution. In the high k limit, the spectral func¬ 
tion dominant contribution is essentially a Gaussian cen¬ 
tered around —A,'^/(2m) and a width {mT)~^^'^k. This 
is explained in Appendix The frequency integration is 
easily done and one finds at this order 


nk 


C'2) 


(16) 


The Tan’s contact is given by 


C(2) = 22 


mV2j'3/2 


,3/2 


r+x, 

[ dxe-^^p2{x) + Q{a-^)Z,„e'^^‘’] 

Jo 


(17) 

where we recognize the contributions coming from the 
contimnun of two-body scattering states and from the 
dimer. It is easy to show that one recovers the known 
result. For instance, at the unitary limit, we find C*21 = 
z^m^T^i/ir [ 21 ]. 

The occupation number Uk is obtained by integrat¬ 
ing on frequency the occupied spectral function A_ (fc, w). 
The results for the BEC regime. Unitary Limit and BCS 
regime are shown in Fig|4(a)| and 4(b) Wc emphasize 





















5 



-200 -too 


(w+//)/ef 


(^) 



-200 -100 0 100 200 

(b) 



-200 -100 0 100 200 


(to+//)/e/- 

(C) 

FIG. 3: The occupied spectral function CfA-( fe, lj) in a 
typical high wavector k limit. T/Tf = 1 and kjkp = 10. 
(a): l/(ifcj:'a) = 1 (BEC regime), (b); 1 /(/cfo) = 0 (Unitary 
Limit), (c): l/(ifcfa) = —1 (BCS regime). The contributions 
with the self-energy of order z is much smaller (blue continu¬ 
ous line) than the contribution of order 2 ^ (red dashed line). 


that it is necessary to include the self-energy up to or¬ 
der in order to get the correct asymptotic beh avior 
rik ~ for large k. We will see in sectioi VI that 

the other contributions of order z'^ do not contribute to 
the contact. 



(^) 



(b) 

FIG. 4: Momentum distribution {k/kpYnk for T = Tf- (a) 
: l/(fcfa) = 1 (red solid line). l/(fcfa) = —1 (blue dashed 
line), (b) : l/(fcfo) = 0. The asy mpt otic values, given by 
Tan’s contact at order (see Eq. jl7| ) are shown as hori¬ 
zontal long dashed lines. For l/ikpo.) = 0, the momentum 
distribution calculated with the term of order in the self¬ 
energy converges to Tan’s contact (red full line). In contrast, 
without the term of order in the self-energy (blue dashed 
line), one does not recover Tan's contact. 


tor”) through the equation 


IV. CONTACT AND PAIR-CORRELATION 
FUNCTIONS INCLUDING 3-BODY 
CORRELATIONS (ORDER :*) 


It turns out that it is more convenient to use the 0- 
pcriodicity of r(P. r). Wc have 


Wc define the pair correlation function 'P(r) = 

— (r, 0) (^'f^'j.) (0,0)). As it was shown in [551 . 

the regular quantity for a short range interaction with a 
bare coupling constant go is P(r) = . 9 o‘P(r). In the limit 
.90 0~, this quantity is related to the two-particle ver¬ 

tex function r(P.r) (sometimes called a "pair propaga- 




(18) 


The contact C is also related to the vertex function |25l 
C = -m^r(r = 0.^-) (19) 


or equivalently, using Eq.i 18 i. C = m^F{r = 0). 
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A. Order z^: two-body correlations 


B. Order : three-body correlations 


In diagrammatic langiiagc, r(P,/3“) is represented in 
Figj^ Aecording to ^ we see dircetly that it contains 
at least a factor corresponding to the global 

time depcndance of pfiir of fields. Indeed, in the Inter¬ 
action representation each annihilation (respectively cre¬ 
ation) operator Ck(T) = cxp(—r(efc — fi)) (respectively 
c|j(t) = oxp(T(eA:— m))) has acxp(pT) (resp. a exp(—^tr)) 
prefactor. We come to the first conclusion that there is no 
such exponential depcndance for close loops, since they 
cancel. The second conclusion is that when a pair of par¬ 
ticle is created at initial time t = 0 and annihilated at 
time T = /3 as is shown in Figl^for the 2-particle vertex, 
we get a factor = ?. Therefore, at lowest 


r 

^^ 

- 1 — 

0 p time 

FIG. 5: Diagrammatic representation of r(P, 


T2 


- 1 -► 

0 P time 


FIG. 6 : The 2-particle vertex r(P, /3 ) at lowest order in the 
fugacity z. 


order, we can tiike the diagrams for F in vacuum, which 
is T2(P.I3~), see Figl^ In the diagrammatic expansion 
this is the sum of idlthe "ladder diagriuns”. We find 
r(P,/3-) = z'^T 2 {P.i.i-) + 0(z^). At order 2 ^, we find, 
using Eqs.lj^ and jlst 




r+oo 

X -f- / dxe~^^p2{x)] (20) 

Jo 


We can perform the integral on P and we get 


At next order, we must find all the diagrams with one 
slashed line. They are shown in Fig 7(a) and 7(b) They 
involve the 3—body problem through T 3 , the 3—body T- 
matrix. Ts can be easily obtained numericafiy, solving 
linear integral cquatioiLS. This is explained in details in 
Appendix B of |5]- These equations have been first foimd 
by Skorniakov and Ter-Martirosian in their seminal work 
on the 3-body problem |2ti| and Rcf.(27|. The analytical 
expressions arc given in Appendix From this we can 
deduce a viriid expansion of the contact 


C = 




(C22^ -f- C 32 ^ + •••) 


( 22 ) 


where the thermal wavelength Aj- = The 

coefficients C 2 (from Eq.( |17| ) and C 3 are plotted as a 
function of the relevant dimensionle.ss parameter Ax/n 
in Fig|^ From Tan’s iidiabatic theorem, we easily find 
that c„ = 3fcn/9(Ax/a), where the b„ are virial ciunu- 
lants entering in the expansion of the equation of state. 
We have cross-hhecked the results obtained for 63 in 
and C 3 (this work) lusing this relation, finding a good 
agreement. 



(a) 



time 


(b) 


i^'^)(r) = (21) 

where is the lowest-ordcr value of the contact given 
in Eq.dl?^. 


FIG. 7: The two diagrams contributing to r*^*(P,^"), and 
hence to C 3 , the virial exp ans ion term of Tan’s contact coeffi¬ 
cient of order 2 ^ (see Eq.J^). 
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A 2 

a 


FIG. 8: The coefficients C 2 (continuous line) and cs (dashed 
line) as defined in Eq.l|22^ for the virial expansion of Thn’s 
contact, in tlie whole BEC-BCS crossover. 


self-energy 

= z'^ F2.i(k,Lj+fi+iO+)+z^ 

(23) 

where 


F2.iik,E) = 


fdP e-^^p2{x) 

(24) 


H2.i[k^E) = - 





E- 


'• 4Tn 


+ 1 - 




2 m 


(25) 


FoUowing the analysis of Appendix|^ we find that this 
diagram does not contribute to the contact. 


V. EXPRESSIONS FOR THE SELF-ENERGY 
INCLUDING 3-BODY CORRELATIONS 


We now show our results for the self-energy at second 
order, i.e. with a double slashed fermionic line or with 
two slashed lines. We give analytical results for < 0 
and equal miisses. We performed niunerieiil calculations 
for the spectriil function at the unitary limit where = 
0 . 


B. 


The diagrams with two are shown in 

Figs 10(a)||10(f) We first detail the calculation of dia- 


grams of Figs 10(a)|lfi(b) 


Figs |10(d)|10(e]| ;md |10(f) 


and 10(c) The diagrams of 


deserve a separate treatment. 
They all involve the 3—body problem through T 3 , the 
3—body T-matrix. 


1 . 


A. 

The only diagram with a double slashed line is shown in 
Figj9| iuid is obtained from FigJ^by changing the slashed 
line into a double slashed line. The analytic expression is 


The diagram of Fig 10(a)| can be redrawn using the un¬ 
folding trick in a way similar to Fig l.')(a) In the time 
domain, we see that it is a convolution product of the 
form (we do not write the wave vector indices for sim¬ 
plicity) 


P-ki 



FIG. 9: The Feynman diagram contributing to E*^* contain¬ 
ing one G*'*’^* line (doubly slashed fermionic line). 

easily obtained from the expression of by multiply¬ 
ing by a factor 2 = and changing into 

which amounts to multiplying the integrand in Eq.(j^ by 
_g-/3rp_k Jiave the contribution to the retarded 


gPrg-(3-r)(e-t-^') / dhdt2e-^'^^T2{h)T^{t2) 

Jv 

X e-"'^'-'^-'=>r2(T-ti - fa) 

where s' = fc'^/(27n), e = (P — k)^/(2m) and V = {ti > 
0, t 2 > 0, r — ti — f 2 > 0}. The double time integral can 
therefore be replaced by the inverse Laplace transform 

^j.r^-(0-r){c+c') f ^e-^*T2(.?-£')T3^(s)r2(,S-£') 

(26) 

where 7 is such that the integrand is analytic for Jf(.s} < 7 
and C-y is a Bromwich contour. In order to get the con¬ 
tribution to the self-energy, wc need to integrate on wave 
vectors, k' is the wave vector of the slashed spin down 
Green's function and P is the wave vector of the two 
T 2 matrices. As usual, wc go to the center of mass 
reference frame, defining the total momentum Pt = 
P -t- k' and relative motion atom-dimer momenta p'j = 
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k' — l/3Pt- In the following, we will need the function 
F{p'i,p'2<s) = t2(s - 3/4(pi)2/m)fJ(p'i,p'2;s)t2(.s- 
3/i{p2)^/m), and the associated spectral function 
/ 03 (p'i>p' 2 ;^) = -1/^^[-F’(p'i.p' 2;^ + iO'*')]- The next 
step is to deform the contour integral in Eq.( |26| iiround 
the real axis, as it is done for the calculation of 
in Eq.(j^. Then, we can take the Fourier transform 
The only dependence on w„ then appears in a 
denominator of the form + fi + e + 6' — x. We can 
safely perform the analyticiil continuation to the real axis 
in order to get the retarded self-energy by replacing iui„ 
by il> + iO'*'. We finally get 

= z'^ F2,2Ak,^ + fi + iO+)+ 

2^-f^2.2,a(A:.w + A( + iO+) (27) 

with a term of order z'^ and a term of order z^. The 
expressions for F 2 , 2 ,a and H 2 , 2 ,a are 


F, 


i,2,^ik,E) = J^dxp^{p\,p\;x) 


g-5e„(P,,p',,k) 

E-x-g^+e.{P<.p\^k) 


(28) 


H' 


■.,2Af^,E} = J ^ dxp3 { p \. p \: x ) 


,-3(*+Sr) 


E-x-^ + e,.{Pi.p\.k) 


(29) 


where 5.(P.p'.,k) = _ 

F 2 , 2 ,a and hf 2 , 2 ,a have the same structure as Fi and 
Hi respectively. The energy denominator corresponds 
to the energy of intermediate states consisting of 3 par¬ 
ticles (relative motion energy x, center of mass kinetic 

energy ^) and of two holes (energy 
(2P,/3-p-,-k)^ 

2 m 


2 . 




The self-energy of diagram Fig l(){b) is denoted 
5 ]( 2 , 2 , 6 ) jjj (.jjp case of equal masses, it is equal to the 


contribution of diagram of Fig 10(a) and we have 


i{2.2.i.) 

-•R 


{k,u) = E',2-2’“>(fc,D.) 


3. 


(30) 


The calculation of the diagram of Fig |10(cJ] goes along 
the same line of reasoning as for Due to the ex¬ 

change of identical fermions, we find a global miniLS sign. 


Moreover, a careful inspection shows that the diagram 
corresponding to the Born approximation for the 3-body 
problem is not a self-energy diagram (not l-particle ir¬ 
reducible) and therefore we should substract it. We de¬ 
note by Ts the new 3-particle T-matrix with the sub¬ 
tracted Born approximation (see Fig |l()(c) ' and p 3 the 
corresponding spectral fimction. For the on-shell 3-body 
T—matrix ts, we have 


^3(Pl-P2;-’) - *3(P'1'P2:S) - ^ 


— s 
(31) 


We have 

E^^’^'"'(A:,cj) =22F2,2,.(A;,u)+A'+f0+)-h2'^H2,2,o(/^'.w+Ai+i0+) 

(32) 

The expressions for F 2 , 2 ,c ^ind h/ 2 . 2 ,c nre 


F'2,2,^(A:,£) = -3^ j^dxp3{p'i.p'2;x 


=.-j3e=(p',,p',,k) 


£ + ^^‘^(P'l-P'2-k) 


(33) 


E2.2A>^-E) = _33y '^P^l^P 2 dxp3{p\.p'2-x) 




£ + ^^‘:(P'l-P'2-k) 


(34) 


where P, = 3(p'i + p '2 + k), £:,,(p'i,p' 2 .k) = 

-f- The factor 3'^ is the .Jaco¬ 

bian of the change of variables. Actuafiy, one can easily 
check that the denominator in Eqs.( |33[ l and ' 34 1 has the 
simpler form E—A:^/(2r7i)—.T-|-((p'i)^-|-(p 2 )^+Pi'P^ 2 )/'^- 


^ j^(2,2,d+i!+/) 

The calculation of the last three diagrams (d),(e) and 
(f) for E(2-2) is more involved. This is explained in Ap¬ 
pendix]^ We also have 

j,^2.2.d+.+/)^^.^^j = Z^F2.2.d+e+f{k.U: + p + iO+) + 


X^H2,2.d+e+f{k,UJ + p + lO'*') (35) 


The results are 


F2,2,d+i: +f{k,E)=-p I 


dP idp'idp '2 


' r +00 /•+“ 


r-l-oo l-l-eo 

I dxi I dx 

Jo Jo 


(27r)9 

g-P(ari + ' * ) P2(xi) 

P2{x- J~~)'^(Pi -P2)(27r)^ + P3(pi-P2;a: 


(36) 
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dPidp\dp'2 

(27r)^ 


^ + 00 /*+3C 

dx 




r + OO r-l 

I (ill / I 

Jo Jo 

P2(xi) 


[^-X-^+Xl + M] (-l-^'l)(-l-^2 

3{P1)^ 


/52(a: - 7 


4 m 


-)^(p'i -P2)(27r)^ + P3(p'i-P2;3;) 


(37) 


whprc e'l 2 are functions of the integration wave vectors. 
Wc have e\ = {p'i)^/m + (Pi)^/(36ni) - • P,/(3m) + 

Pj - k/(m) — k • P(/(6m) + k^/(^m) and an equivalent 
expression for e'j. 


P-kl P-kl 



0 r' t" P time 0 t' r" 0 time 

(a) (b) 



FIG. 10: The 6 diagrams contributing to — t'), 

including the 3-body problem. 


VI. NUMERICAL RESULTS IN THE UNITARY 
LIMIT 

Our numerical calculations were performed in the uni¬ 
tary limit. This brings simplifications, since one can 
rescale the wave vectors entering in ig, and there arc 
only two variables in tj, which significantly decreases the 
amount of computation. It is one of the reasons why 
we studied numerically the unitary limit. The imaginary 
part of the self-energy is calculated directly numerically, 
and the real part is obtained using Kramer.s-Kronig rela¬ 
tions. 


A. High negative frequency, high wave vector 
behavior and Tan’s Contact 


We first consider the large momentum k and large neg¬ 
ative ui limit of the self energy, where we expect to recover 
the physics of Tan's contact. In Fig ]Il| is shown the imag- 
iniuy part of the self-energy at order (the imaginiuy 
part of the sum of H 2.1 and hf2,2,•••)• In agreement with 
the lowest order calculation, we find that ^2,1 find ^2,2, 
do not contribute. The niunericid results show that the 
main contribution comes from H 2 . 2 ,a ^nil H 2 , 2 .b- Mid can 
be well fitted by a Gaussian. Moreover, by integration on 
oj, we get the occupation number n^, from which we can 
find a value for the virial expansion cocfiicient C3 (see 
Eq.i22i) of Tan's contact at order z^. Wc find a good 


agreement (within a few percent) with the value foimd 
in section |IV B| at unitarity, using a virial expansion of 
the two-particle vertex F. We expect this featme to be 
valid at any order of the expansion in fugacity z. 



FIG. 11: The imaginary part of the self-energy at order 
(see text) for fcAr = 20\/7r- 


B. Spectral function A{k.uj) : effect of three-body 
correlations 

We come now to the main results of this work on the 
effect of three-body correlations on the spectral function. 
Wc have calculated numerically the retarded self-energies 
corresponding to diagrams of Figsj!)|and |l()[ In Fig |12(a') 
the imaginary pmt of (^2,2,0 the contribution o' 

order is plotted for kAr = 2^/^ and 2 = 0.3. Using 
the virial expansion of Appendix [a| up to order 2^, we 
find T/Tp = 1.60. Wc observe a sharp structure at the 
point uj + fi = k'^/{2m). We have found a sharp varia¬ 
tion of the form ■\/\ui -f p — A:^/(2m)|. This is explained 
in Appendix Wc foimd a more pronounced variation 
at w -h p > k'/{2m) than at u -|- p < k^/{2m) because 
the cocfiicient multiplying the square root is larger in the 





































































10 


former case. Indeed, for the ehaonel I = 0 angular mo¬ 
mentum ehannel, this ratio is about 8. This singularity 
is associated to the 2-body matrix t^is) and as a result 
this singular behavior is observed in other ternrs of the 
second order self-energy. In Fig 12(b)[ we plot the imag¬ 
inary part of the total self energy of order 2^ for the 
same wavevector and fugacity (including 0(2^) contribu¬ 
tions of diagnims of Figs and 10 i. Here too, we can 
observe a sharp structure aro\md w + = k^/{2m). We 

also observe that the contribution to the imaginary part 
(see Fig 12(b) ■ of the self-energy is positive, whereas the 
self-energy must have a negative imaginary part. Actu¬ 
ally, this is consistent with Fig |ll[ where the imaginary 
part of the self energy is also positive, consistent with 
the fact that cg is negative. However, this restricts the 
rciusonablc values of the fugacity since, for example, for 
2 = 0.5 and fcAx = the spectral function can have 

some negative contributions. 

The second order spectral fimction (including all 0{z^) 
terms in the self-energy) is plotted in Fig |13| for 2 = 0.3. 
Given the features of the imaginary part or the .second 
order self-energy, we find a very different behavior com¬ 
pared to the lowest order calculation 1121 . Similar to low¬ 
est order, there exists a two-peak struetme at small mo¬ 
mentum, which gradually evolves into one peak at large 
momentum. However, the two peaks are very asymmetric 


and the left peak dominates. This is shown in Fig 13(a) 
for 2 = 0.3 and kA^ = We conclude that 3-body 

correlations can be quantitatively important for the spec¬ 
tral function, already at a high temperature correspond¬ 
ing to 2 = 0.3. Naturally, we expect this will become 
more pronounced when the temperature decreases, and 
the fugacity increases. Our calculations imply that it is 
not sufficient to include only 2-body correlations in order 
to obtain the correct spectral function, in piuticular close 
to the critical temperature. 

Concerning the issue of convergence of our results with 
respect to higher powers of the fugacity 2, the only thing 
we can safely say is that at this value of z (= 0.3), the 
virial expansion of the equation of state at third order re¬ 
produces accurately the experimental results of Ref.[2]. 
Hence, we may expect that the expansion has also con¬ 
verged for the spectral function. 


VII. CONCLUSION 

In this work, we have extended a high temperature, 
virial-like. expansion to a dynamical quantity such as 
the self-energy or the spectral function. Even though the 
method can be applied to bosons as well as fermions, we 
have considered in detail the BEC-BCS CTOssover for spin 
1/2 fermions in 3£5. Wc have calculated the self-energy 
on the real axis up to order 2^ (2 the fugacity), includ¬ 
ing thret'-body correlations, for < 0 (Unitary-BCS 
side). We have performed numerical calculations in or¬ 
der to evaluate the spectral fimction in the imitary limit, 
including three-body correlations. We find a quantitative 



(ci>+p)/T 

(^) 



(a)-hp)/T 

(b) 


FIG. 12: (a) Imaginary part of w) for fcAj- = 2^ 

and z = 0.3. (b) Imaginary part of Sl^*(A:,w), the second 
order self-energy, for kAr = 2</Tr and 2 = 0.3. 


importance of three-body correlations even in this high 
temperature regime. We have also calculated the third 
order (0{z^)) virial expansion of Tan’s contact, in the 
whole BEC-BCS crosso^r. 

The problem of the mixtme of fermions of different 
masses is a natural extension of our work. Experimen¬ 
tally and theoretically, this hius been studied in |28) with a 
mixture of ^Li and atoms. There, we expect 3-body 
correlations to become particularly important since, at 
least on BEC side, one expects |2iJI a three-body boimd 
state in / = 1 (iiannel for mass ratio around 8. 
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(ci)+|i)/T 

(^) 


At fa from the two other dimcnsioolcss parameters T/Tp 
and l/{A.'fa). 

The Fermi momentum kp is defined siieh that the total 
density n = (^.•^)/(37r^) and the Fermi temperature is 
defined by Tp = fe|./(2m). Ax = ^27r/(m T) is the 
thermal wavelength. The virial expansion up to order 
states that nA^ = 2[z + 262(Ar/a)2^]- From this, wc 
get the two equations 


At 

a 


Z + 2b2{^)z'^ 

a 




s-v/tt Vr^/ 


(Al) 

(A2) 



((D+h)/T 

(b) 



FIG. 13: Spectral function at second order (blue solid line), 
compared to the lowest order ones (red dashed line).The fu- 
gacity is 2 = 0.3 [T/Tp = 1.60). (a) kAr = v^/5, (b) 
kAr = (c) feAr = 4y^ 


For given values of l/(A:f a) and T/Tp, we find ^ from 
Eq.( |Al[ ). This gives the virial cumulant fe2 (Ar /Q) (see for 
instance [S| for its expression (3) in Eq.(A2i, which for 
given TjTp becomes a (second order) equation for the 
fugacity z. In this way, if TjTp = 1, we find s = 0.17 
for l/{kpa) = 1 (BEC side), z = 0.5 for l/{kpa) = 0 
(unitary limit) and ? = 0.7 for l/{kp a) = — 1 (BCS side). 


Appendix B: BCS limit 

The BCS limit is a weak coupling limit, therefore 
we expect that the real part of the self-energy will 
be given by a mean-field like term gn^, where g = 
ipajm is the coupling constant and at lowest order 
ni = 2(mr/(27r))3/2. This yiek is an energy shift 

y^(2m)/7ra z Noti ce that this shift is negative, as 

can be seen in Fig |2(e)[ Wc can recover this results us¬ 
ing the explicit expression of the self-energy Ti{k,U)) in 
the limit a~^ —> —oo. The width of the peak is ob¬ 
tained using a Fermi Golden Rule approach. The ini¬ 
tial state is obtained by the creation of a spin '[ fermion 
in a gas of 2N classical fermions (N spin N spin 
;). Let us denote the momenta of the spin 4- fermions 
Pii • • ■ - Pn- The created fermion scatters on a spin 4- 
fermion and the final states are obtained by the cre¬ 
ated fermion with a momentum k -|- q while the scat¬ 
tered spin 4- fermion, for instance the one of initial mo¬ 
mentum p/y has now a momentum pjv — q, by momen¬ 
tum conservation. All the other particles keep their 
momenta. The difference in energy in this process is 
(k-t-q)^/(2m)-f(pAr-q)^/(2m)-(A:V(2m)-fp^/(2m)) = 
g^/j7i-|- (k — pjv) • q/m. The inverse of the lifetime of the 
initial state is therefore 


Appendix A: Finding the fugacity for given TjTp 

and \l{kpa) 


^ = 2 ^ 9 ^ j (g^/m-(-(k - pw) ■ q/m)) 


1 

47r 


5 ^ 772 |k - p/v| 


(Bl) 


The goal of this section is to explain how one can find 
the two dimensionless parameters z (the fugacity) and 


The next step is to average on the fermions state. iLsing 
a Boltzmann distribution. The calculation is similar to 
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the one of Appendix |(' 1[ and we just give the residt 

1 Iza^ '•+^ 


Tk 


TT 0k 


L 


dpp^ 


. {K-py <■ (>+?)■= 

’2m _ e 2m 


;b2) 


This result can actually be recovered from the expression 
of the imaginary part of £(A;, w) at lowest order in z in 
the limit a~^ —^ —oo, (or u + p = k'^/{2m). Indeed at 
lowest order in the fugacity z, the imaginary part of the 
self energy is the imaginary part of w + ^ + iO"'"). 

In the BCS limit a~^ —oo, and we can replace P 2 {x) 

in Eq.(j^ by In a perturbative approach, 

we can aLso replace ^ by k^/(2m) and wc find for 
—Q(E(fc, k^/{2jn) + iO'*')) at lowest order in z and a 


rzl- --^e P 2m / 

y(27r)3 Jo 


+ 00 
dx- 


ll/2' 


‘^S{^-X+ 
2m 


{P - k)2 p2 


2m 


4m 


The integral on x can be performed, and by compari¬ 
son with Eq.| |B2[ f. we find the expected result (after the 
change of variable P' = 2k — P) —23(S(A;, fe^/(2m) + 
i0+)) = 1/Tk. 

When wc compare our niunericid results with the per¬ 
turbative results obtained in the BCS limit, we sec that 
we need to go quite far in the BCS regime in order to get 
a good quantitative agreement. 


Appendix C: BEC limit 

1. Occupied spectral function of a gas of classical 
dimers 

We recall the definition of the occupied spectral func¬ 
tion A_(fc,a;) 


A_(fc,w) = ^^"|(m|ck,tl'2)l^^(^ + Em - 

n,m 

(Cl) 

where the sum extends to all the eigenstates |n) and |m) 
o( H' = H — pN where H is the Hamiltonian of the sys¬ 
tem and N the particle number operator. is the grand 
canonical partition function. In the BEC limit, the eigen¬ 
states |n) are given by t-L dimers in states of momenta 
Pi, • • • .Pa/ (states with an odd number of atoms have 
much higher energy since they don’t benefit of the bind¬ 
ing energy). For a given dimer niunber M, the action of 
the annihilation operator Ck,f couples to states consi.sting 
of a single 4- particle of momentum Pj — k (i = 1, - • • M) 
and dimers in states Pi ,•• • , P,_i, P,+i. • - • Pa/. There¬ 
fore we have E,n - E„ = p + {P; - k)^/(2m) + Et, - 
P)^/(4m). The indistinguiahabOity of the dimers makes 
that there is a factor M (the niunber of possible val¬ 
ues of i) in the summation. We choose Pa/ the mo¬ 
mentum variable of the dissociated dimer. The matrix 
clement in the dilute limit is (we basically ignore the 
effect of the — 1 dimers) (m|ck,t|n) — (Pm — k j. 


|rk,-|.|£l : Pa/), where the dimer state of momentum 
Pa/ is \D : Pa/} = Eq V’q4„/2+qt4„/2-qtl0) (|0) iH 
the vacuum state), ipq = (87r/a)^^^l/(g^ -|- a“^) is the 
dimer wave fimction. We find for the matrix element 
(m|ck,-f|n) ~ ‘(^k-PAi/2' The summation on the AI — 1 
dimer variables give the partition function of A/ — 1 cliia- 
sical particles of mass 2m. It factorizes from the rest 
and is X];a/>i M/M\(J2p e-^(^'/(4m)-E,-2p)jA/-r ^ 
which simplifies with the at the denominator. Wc 
find 


Stt i^(t^ -f + Es - ^) 

a (a-2-f (k-P/2)2)^ 


This can be written, due to the Dirac fimction 


Stt 


A_(fe.w) = 2^—^e'- 
m^a 


1 




/( 


dP 


e ^S{uj + p + 


^ + ' 2f: 

(P-k)2 , 




(2?/)^ ' ' ' 2771 

From this expression, we easily find the tluesliold 
U + p < —Et + fc^/(2777) 


The integral in Eq.(C2i can be performed analytically, 
and we find 


A_{k.Lj) 


gpl _apl 
^2 e im -e ^ 

Trak + 


(C3) 


i( u + p < —Et + k^/{2m) and 0 else. Wc have 
P_ = 2{k — y/?n|u; + p + Et — fc^/(2m)|) and P+ = 
2{k + + p + Et - k^/{2m)\). 


2. Particle peak 


At lowest order, the situation is similiu to the BCS 
limit. One finds a positive mean field like shift of the 
peak located around lo + p = k^/{2m). and a broadening 
given by the Fermi Golden rule like calculation Eq.i B2 i. 


Appendix D: Limit of u —> —oo, fc —» oo 

We calculate the dominant contribution of the imagi¬ 
nary part of — I/ttEi (A:, w + ^ -|- i O'*") in the limit of high 
k and large non pasitivc frequency w -f p. Due to the 
Boltzmann weights and in this limit, 

we find that P ~ (777/^)1/2 X < {ml^y^. This 
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means that in this limit -4^10) + and 

I -C |u) + fi\, /m. Hence we find 


- '^\Hi{k,LJ + fi + i 

TT 


C'S{lj + /.i 


.0^)1. / 


dP 


(27r 


3 

P k. 


'4m X 


2771 


m 


where the constiint C' is given by 

r+oG 


C = f dxe-^V2(x) + e(a-^)e^^‘'Z„ 

do 


(Dl) 


(D2) 


The reason why we do not neglcet is that k is large 
in this limit Jind this sejihir product can be comparable 
to w + ^ 4- k'^/{2m) in the relevant case where w + ~ 

—k^/{2m). The integral on P can be performed exactly 
and we find 








0 T. 


0 T, 


0 T, 


0 r. 


^2 

(a) 


(b) 


/3 


time 




time 


time 


fi ti 


time 


which is a Gaussian centered around w + /r = —k^/ (27n) 
and of width of order {T/m)^^^k. 


(D3) 

FIG. 14: Unfolding trick : the slashed line — ri) can 

be drawn as the product of two unslaslied lines, (a): 

T 2 > n, (b): T 2 < n- 


Appendix E: Unfolding trick 

We show here how a slashed propagator r) 

can be reinterpreted in terms of Feynman diagrams. We 
consider a general shxshed fermionic line between two 
times Ti and T2 with 0 < Ti < ^ and 0 < Ta < ^. We 
havcG'°-^'(p,T2-ri) = = G<°'°’(p,^- 

ri)G'°’°*(p, T2)- Therefore, a shxshed line connecting the 
time Ti to the time T2 can be replaced by two imslashed 
G<o,o) lines. The first connects the times Tj and /3 while 
the second line connects the times 0 and r2. This 

is shown in Fig 14(a) and|14(l7) 


Appendix F: Calculation of r*^*(P,/3 ) 


The diagrams of Figs |7(a)|7(b) can be drawn iLsing 
the ’’unfolding" trick (see App<‘iKiix as shown in 
Fig.s{^ W (b)| The analytical expression in the imagi¬ 
nary timeoonmn involves the intermediate times t\ and 
T 2 iis shown in Fig{Tg Ti is the time of the last inter¬ 
action in the incoming T 2 . while T2 is the time of the first 
interaction in the outgoing r2. The whole diagram, e.g. 
Fig |15(a')] is ii convolution product whicii ciin be written 
schematically, without wave vector indices, 

[ dUdt2e-^'^'T2{h)n{t2)e-^^^^-‘'-*^^T2{0-h~t2) 

Jv 

where ti = ti and fg = T2 — ti and the integration do- 
miiin 2? = {(fi.fg) : ti > 0.t2 > 0,^ — ti ~ t 2 > 0}. 



0 T, 


0 X, 


X, |5 


(a) 



(b) 


time 


Xt (5 


b'me 


FIG. 15: The two diagrams of Figs 7(a)|7(b) drawn in the 
unfolded way. 


This convolution product can be written iLsing Laplace 
transforms |^. For the diagram of Fig 1.5(a) we find 


xT3,t[(k,£fc),(k,£fc);(k-hP..9)] 












































14 


III the Bromwich contour C-y, 7 is such that the integrand 
is analytical for 3^(.s) < 7. The summation on k is done 
by going to the center of mass reference frame of the 3- 
body problem as is explained in |^. Ta is the 3-particlc 
T-matrix as defined in terms of Feynman diagrams in 
[57]. We find for 


X 



9 mpi 


3 fl(i.i+p>^ 

e g - - (pi -> -Pi) 



Appendix G: Calculation of £(2,2,ti+c+/i 


The summation of the three diagrams of 


Fig s|10(d)|l(l(e] and 10(f) is shown in Fig 16(a 


to 

Fig |16(d)| TK? ne we see that the sum of the three 
diagrams can be factorized. The first term in the first 
parenthesis of Fig 16(d) can be simplified. Indeed it 
contains the product of two T 2 imd two If wc 

denote by T2 the time at the exit of the r2 on the left 
and Ti the time at the entrance of the T 2 on the right, 
we have the contribution 


l'{P;T2,'ri) = J j dxidx2T2iP.xi~Ti) 

xG<°'^>(k.i2-xi)G‘°'^>(P-k ,l2-Xi)r2(P,T2-X2) 

(Gl) 


three terms, and the integral must be down on the imion 
of the three domains T> = Vi \JT >'2 Wc easily find 

V = {(fl,f2),^l > 0,f2 > O,0-T-ti-t2 > 0}. Therefore 
the sum is a convolution product and is equal to 


/;27rf 4m (s-ei)(s-£2) 


(G4) 


Notice that it is a fimction of the time difference r’’ — r'. 
as it should, but it is not obvious. In order to get the 
Fourier transform, wc deform the contoiu along the real 
axis. By doing so, we must take into accoimt of the 
branch cut of t 2 {s] and the two poles at s = £1 and £2. 
For a~^ < 0 wc find 


r+o® g-(3-T)ii p2 

- £i)(ii-£2)'°^^^^ 4m^ 


'/ 


+ 


47r f e-l/S-i-lEi N 


where V is the Cauchy principal part. Notice that at 
imitarity. the second term vanishes. 

Inside the interval [r', r"], wc obtain the two terms in 
the second parenthesis of Fig 16(d) The first term, for 
given wave vectors P. ki and k2, is equal to 

-T2(P-ki +k,T)G(°'°>(k,.r)J(ki -k2)(27r)3 = 
e-^-u^r2(P - ki + k, r)<S{ki - k2)(27r)^ (G6) 


This can be written in term of inverse Laplace transform 


where V = {(xi,i2)-0 < X2 < T2 < rj < xi < 0}. It is 
shown in Appendix that 

/{P;T2,Ti) = -r2(P,^ + r2-n) (G2) 


The sum of the three terms in the first parenthesis of 
Fig 16(d)|Ls cqiml to (for given wave vectors ki and k2) 


f ^e-^^T2{P - ki + k,.s - £fc.)<S(ki - k2)(27r)3 

(G7) 

The second term is of the form (we do not write all wave 
vectors for simplicity) 


f dtidt2r2(P,/3-r-fi 

Jvi 

JVi 

- f dtidt2r2(P.t2)e-'^‘‘i+®=f^-'"-‘i-'^» (G3) 
Jv-, 


f dtidf2e-^^i'^Ta(fi)ra(t2)e-^^^<^-‘‘-‘')r2(r-fi-f2) 

(GS) 

which is again a convolution product and can be writ¬ 
ten as the inverse Laplace transform of the product of 
Laplace transforms 


where £1,2 = (fci,2 + (P~ki,2)^)/(2m), r = t"-t' and *1,2 
denotes time differences. In the second term, we make 
the change of variables = 0 — t — ti — t 2 and = ti. 
Wc get a new integration domain 2?2- third term, 

we make the change of variables t '2 = 0 — r — t\ — t 2 
and fj = ti, and wc also get a new integration domain 
'D'^. By doing so, wc find the same integrand for the 


L 


ds 

2ni 


e-^^T2{s - £fc,)T3(s)r2(s - £fc,) 


(G9) 


The next step is to deform the integration contours of the 


variable s in Eqs.( G7 ■ and (G9 ■. This is done as before 


by introducing the spectral c 


cnsitics p2 and p.^. The final 
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result for the sum of the diagram (d)+(e)+(f) is 

= -z’^V j 


P-k.T 




P2{x'i) 


(i'i-ei)(xi-£^) 





U^icEeMEi^rsEi 


k^r bjt 


P2{x' ~ - P2)(27r)^ + p3(p'i,P2:a;') 

4 771 

(GIO) 

where we have made the changes of variables i' = i — 
P^/{Gm), x'l = xi -F^/(477 i) and e'j 2 = £1,2 --P^/(4r7i). 
If we take the Foiirier triuisform and make the analytic 


continuation, we find Eqs.(|36[), (37 . 


Appendix H: Proof of Eq.|]G2 


The diagram of Fig 17 is denoted /{P:r”,r'), with 
r” < t' . The two-particle vertex r(P.T) has an explicit 
imagiiiiiry time dependance due to the time evolu¬ 

tion of the two created fermions. We define 


r{P,T) =e-2^"r(P,r) 


(HI) 


This new vertex can be expanded in terms of the fugacity 
r(P,T) = En>o-2''f*'‘>(P,T), where 2 = Moreover 
we have the periodic boundary condition r{P,r < 0) = 
r(P,/3-h r). which yields f(P,r < 0) = 2^f{P,/3-H r). 
Hence we have f<’‘+2'{P,r < 0) = + r). For 

n = 0. it Ls easy to see that ff°>(P.r) = r2(P,T), the 
2-particles T-matrix in vacuum. Moreover, it is also easy 
to show that for r < 0. -f<2){P,T) is precisely given 
by the diagram of Fig |17| (for r < 0, one needs at least 
two slashed fermionic lines which go backward in time). 
Therefore we have 



IP-k, P-k, 

(krElkpEfc 

+ 

p-k, P-k, p-k,l 


0 t' r” 

i time 

y 0 time 


(b) 



0 r' r"Sii"e 


kt +. 




P-ki p-k; 

= 


;P-k2P-k. P-ki! 


0 r' r’'iliin>« 0 

(C) 


+ i'«2 ^ k _ '• ^ 

I **1 ' kj r 5 
-J-^—u 

j' t" pVn* 



0 r* r" tf ^ ^ 

P-k,T 






kii 

kji 

+ 

>7k: 

“k(T 

iP-k, 

ikiriTt 




P-k,i 

0 r 

I 

' 0 time 

t' r" 


T 

' ? 


/(P; r , r') = -T2{P, ^ + r" - r') 


(H2) 


p-k, 


• f k, ri-s • "(k] - k2) + k; r H— J 


P p — kj 


0 1' 1" 


0 t’ r'' ^tJme 

(4) 


Appendix I: Structures in 3(E*^'^'“*(fe,Lj)} 


In Eg. | |2S[ ). we can make the change of variable Pj —> 
P' with P, = 3P'-t-3k-|-3/2p'i. In this way the imfiginaiy- 
part of S^^'^ “*(A;,(j) is given by 


FIG. 16: Diagrams of Figs infdllKlf^ and |10(f)| written in 
the unfolded way. (a) is for E^’^' . (b) is for E^^ and (c) is 
for E^’^’7. (d) shows how they can be factorized. 


pi 2 g \2 

X p3(p'i,p'i;x)«5(w-x+ --I- - 7 ^^) (II) 

m 4 771 

with w = u)+ii—k^/[2m). The integral on x is performed. 
Due to the Dirac function, the integral on x can be per¬ 
formed. We find 1 = 0; + 3/i{p'i)^/m + P'"^ jra = xq > 0. 


In this way, we find 

- i (t,,.)) = 3 ’.= I m (P',. P'.; x„) 

P'2 3 (p', )2 

x0(w-t--+ 7^^) (12) 

m 4 m 
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P-k 3 T 



Ill the iinitfiry limit, p3(p'i.p'j:xo) is given by 


^’aCp'i-p'ilxo) = —7p7Tr3'[i3(p'i.p'i;a:o)] 

rM /. I j_ \ _ I 


m‘- j:,+ 

Z)I\2 


_;r5(w+^)K[t3(p'i.p'i:xo)]} (13) 
in 


The first term (respectively the second term) gives 
a contribution to —l/7rQ(E*^'^’“*(fc,u;)) denoted 
/j(ij,fc) (respectively Hence we have 

-l/7r3(S(2h,a)(fc^^)) = k). 

We first consider l 2 {io,k). Due to the Dirac function 
(5{(I> + ), we see that this term vanishes for w > 0. 

Moreover, for w < 0, the integration on the norm of P' 
gives a term, for small |tj|, of order dP'{P')'^5{—\Lj\ + 

~ y/WI- come to the coneliLsion 

/2(d;,A;) = 0(-£lj)-\/|ih[j^(£l». fc) (14) 


where P is a regular function of Oj and k. This residt is 
of course confirmed by our numerical calculations. 

In the expression for Ii{u>,k), w enters in the imagi¬ 
nary part of ts. in the 0 fimetion in Eq.(jl^ and in the 
prefactor l/(u)-|-P'^/m) in Eq.(j^. We first consider the 
uj dependence from this last term. We have found argu¬ 
ments which show that the other w terms will be of order 
wln(£i;) and lo, which are higher order terms. Following 
this method, we find for the difference /j (w, k) — /i(0, k) 
(we take the mass m = 1) 




■ + {P'y 


X [*3(p'i. p'l; 3/4(p;)= + P'^ + i^)] (15) 


We see from this expression that for small Cj, the inte¬ 
gration on P' will be limited to values of order 
At lowest order, we can set P' = 0 in the expression 
in factor of l/(£l’ -f (P')^)- The integration on P' gives 
V ffj dP'l/(ii; + P'^) = 7r/2\/59(w). As a conclusion, we 
find 


7i(£l;, k) — 7i(0, fe) = -1- (16) 

The prefactor C (which depends on k and ^) can be easily 
calculated once ts is known. The next order terms of 
order il/ln(tl)) and ui are fitted to reproduce the numerics 
with reasonable cocfiicients of order unity. We find that 
all the arguments of this Appendix give good analytical 
supports to onr numerics. 
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